Partially disordered state near ferromagnetic transition in MnSi 
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The polarized neutron scattering in helimagnetic MnSi at low T reveals existence of a partially 
disordered chiral state at ambient pressure in the magnetic field applied along (111) axis below 
the first order transition to the non-chiral ferromagnetic state. This unexpected phenomenon is 
explained by the analysis of the spin-wave spectrum. We demonstrate that the square of the spin- 
wave gap becomes negative under magnetic field applied along (111) and (110) but not along the 
(100) direction. It is a result of competition between the spin-wave interaction and cubic anisotropy. 
This negative sign means an instability of the spin wave spectrum for the helix and leads to a 
destruction of the helical order, giving rise to the partially disordered state below the first order 
ferromagnetic transition. 

PACS numbers: 75.25.-|-z,61.12.Ex 

Non-centrosymmetric cubic helimagnets such as MnSi, FeGe, FeCoSi are the subject of the intensive experi- 
mental and theoretical studies for the last several decades. Their single-handed helical structure was explained by 
Dzyaloshinskiii. The full set of interactions responsible for observed helical structure (Bak-Jensen model) was es- 
tablished later in^ii^ in agreement with existing experimental data (see for example^ and references therein). The 
renascence in this field began with a discovery of the quantum phase transition to a disordered (partially ordered) 
state in M nSi at high pressure- and- . The following properties of this state attract the main attention: i) non- Fermi- 
liquid conductivity, ii) spherical neutron scattering surface with the weak maxima along the (110) axes^-f^, whereas 
at ambient pressure Bragg reflections are observed along (111)^. These features and the structure of the partially 
ordered state were discussed in several theoretical papers (see^ii^iiiii^ and references therein). It should be noted also 
that the spherical scattering surface with maxima along (111) was observed at ambient pressure just above critical 
temperature Tc ~ 29K. This experiment was explained using the Bak-Jensen modeli^. 

These studies shadowed an important problem of the helix structure evolution in the external magnetic field H. In 
particular, simple phenomenological^'* and microscopical^''' theories predict the smooth second order transition from 
the conical to the ferromagnetic state. The spin component of the cone parallel to the applied field is proportional 
to the magnetization and it increases as {H/Hc) up to its saturated value. This prediction is in agreement with 
experiment''. The perpendicular, rotating spin components fade away with the field and according to the plain theory 
the helical Bragg reflections must decrease as {H^ — H^)/H^ where He is the critical field for the ferromagnetic 
transition. This, however, contradicts to the experimental facts if H || [111] (Seei^ and Fig.l). The experiment shows 
that the transition is of the first order and the Bragg intensity does not follow the law given above. 

In this paper we demonstrate that although from a general expression for the ground state energy one can expect 
the second order transition at the critical field He, this is true for H || [100] only but the situation changes if H || [111] 
or H II [110]. In the last cases the spin-wave spectrum is unstable in the field range Hi < H < He due to the cubic 
anisotropy as the square of the spin-wave gap becomes negative and the helical long-range order decays in this field 
interval. The ferromagnetic state occurs above He ■ Hence we have a region below He where the partially disordered 
state coexists with the almost saturated magnetization. It is a reason why this state has not been noticed in the 
earlier macroscopic measurements^,^. 

Let us consider the conical helix with the lattice spin Sr = 5'rCr + •S'r^^r + S'rCr where 

Cr — c sin a + (Ae'*' -^ + c.h.) cos a; 

m = z(Ae'''-^ - c.h.); (1) 
Cr = c cos a — (Ae*'''^ -I- c.h.) sin a, 

where A = {a—ib)/2 and the unit vectors C, r) and ^ form the right-handed frame. If a = we have a plain heli»i^. The 
spin operators are given by the well known expressions: S^ = S— (a+o)R, S^ = — i(S'/2)-^/^[aR — a^— (a+o^)R/(2S')] 
and si = (5/2)i/2[aR + a+- {a+a^)R/i2S)]. 
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Similar to inJ^, we use the following Hamiltonian 

+ I] Sq.iS-q.iqf} + K J2 Sci,iSp,iSh,iSf,i (2) 

l—x,y,z l—x,y,z 

where the first term is the ferromagnetic exchange interaction, the second term is the Dzyaloshinskii interaction (DI) 
responsible for the helix structure^. The following terms are the anisotropic exchange, the cubic anisotropy and the 
Zeeman energy. They determine the orientation of the helix vector k in the magnetic field (seei^r^ andii). The 
following hierarchy of interactions holds^i^: Jq >> Dg/a >> Fq/o^ ^ K where a is the lattice constant. Replacing 
Sr — > S'Cr one gets the classical energyi^ 

Fri = b { —— H 
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-SDoCk • [a X cos^a + SH^^ sin a + Ecub 



where A = S{Jo — Jkj/k'^ is the spin-wave stiffness at q >> k, L(k) = ^ kf{af + bf) and is the field component 
parallel to k. 

Using (3) and (4) and taking into account that the single ion contribution Ecub does not depend on k we get^^ 

k = kiDo/\Do\)[ci - SFoiaf + bf)/{2A)] 

Eci = -{SAk^/2)[l-SFoL{c)/{2A)]cos^a + Hnsma + Ecub, 

where k — S\Do/A\ and L{c) is a cubic invariant. For Dq > or Dq < we have the right or the left helix, 
respectivelj^. For Fq < (>)0 the helix vector k is oriented along the (111) ((100)) axes as the invariant L has two 
extrema 2/3(0)^. 

Neglecting Ecub in Eq.(4) one obtains the conical state with sin a = —H\\/Hc if H\\ < He where He = Ak'^(l — 
SFoL/{2A)] ~ Ak'^ and the ferromagnetic state for H\\ > Ho-^^. It can be shown that Ecub gives a negligible 

contribution to He and to sin a. The principal parameters of the magnetic structure for MnSi are: A ~ 52meVA^ , 
k ~ 0.038A and He ^ 0.6T ~ Ak'^ in agreement with the theory (seei^). Another energy Fofc^ O.OlmeV" ~ O.IT 
was estimated from the anisotropy of the critical neutron scatteringi^ and from the reorientation of the helix axis in 
the magnetic fieldi^. The value of K will be estimated below. 

The classical energy depends on only. The general expression for the ground state energy, derived in^'''. It 
contains yet another term of purely quantum nature, which depends on the spin- wave gap A as a parameter. As shown 
ini^ if H± > A^/2 the helix wave vector k is directed along the magnetic field. For MnSi A ~ 12 fj,eV ~ O.IT ^ He^- 

In a linear spin-wave theory the gap appears due to the cubic anisotropy but it equals to zero at — 02£. There 
is yet another contribution to the gap, which is a result of the spin-wave interaction. We begin with the former. For 
evaluation of A^ one has to consider the uniform part of the bilinear Hamiltonian, which is given by 

Ho = Eoa+ao + Bo{al + a+^)/2, (5) 

and Aq = Eq — B^. If one neglects the cubic anisotropy at iJy < He then one has Eq ^ Bq ^ {He/2) cos^ a and the 
gap is zeroi^fi^. Taking into account Hcub, one obtains after simple but rather tedious calculations 

Eo = [He cos^ a)/2 + Ei + E2; Bq = (He cos^ a)/2 - E2 
El = -4A[(1 - L) sin* a + 3L sin^ a cos^ a 

-I- (3/8) (2- L) cos" a]; (6) 
E2 = (3A/4)[4Lsin2 a + {2 - L) cos^ a]; 
^I'ub = {Hecos'a + Ei){Ei + 2E2), 

where A = S^K. These expressions hold at H\\ < He and, indeed, A^^^ = in zero field only2£. 

For H > He we have cos a = and if without cubic anisotropy Eq = H — He , Bq = and the gap A = H — Hcr^, 
while if taking it into account the cubic anisotropy we get 

Akb = [H~He- 4A(1 - L)][H - He + A(10i - 4)]. (7) 
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For the two principal directions Lm = 2/3 and iooi = and we obtain: 

^{H-Hc- 4A/3)(F -Hc + 8A/3); 



^Ctif), [1,1,1] 
^Cu6, [1,0,0] 



Akui.o.ol ^{H-Hc^ 4A)2, 



(8) 



Thus, one comes to the important conclusion: A^^^^ ^ ^-^ is negative ai H = He for both signs of K^^ . This 

circumstance is decisive for the stability of the system if one takes into account that the contribution to A^, stem 
from the spin- wave interaction, is proportional to cos'* a and disappears at H close to He- 

Let's consider the spin wave interaction to the gap Aj„j. As the DI breaks the total spin conservation law it must 
lead to the spin- wave gap. However, in cubic crystals this interaction is very soft [see Eq.(2)] and the gap appears as 
a result of the spin-wave interaction only similar to the case of pseudo-dipolar interaction in antiferromagnets^^)^. 
In the one-loop approximation it consists of both the Hartree-Fock (HF) part evaluated ini^ at -ff = and the 
second order contribution from three-point spin-wave interaction, which appears due to the helical structure. It was 
ignored in^^. The diagrams for both contributions are shown in Fig. 2 where lines correspond to Green functions 
Gq = — < Tflq, a+ > and i^q = = — < Toq, a_q >, which in tj-representation are given byi^: 

= + ; F^{iuj) = - , . ,f 3 , (9) 

where = S{Mo ~ Afq) + Bq; Bq = (5/2)(Mq - Jq) cos^ a; A/q = J+ + 2Dq(k • c), J± = ( Jq+k ± Jq-k)/2 and the 
spm-wave energy eq = (E\ - B^y/'^. Although at q < 1/a these expressions give the same result as obtained iiiii: 
Eq = Aq^ + Bq] Bq = {Ak^ / 2) cos^ a and eq — Aq{k-^ cos^ a + q^Y^"^- However for the present consideration we 
need them for all q as the formulae for A|,^j contains the sums, which saturate at g ~ 1/a. 

The contribution of the forth-point interaction to A^ was analyzed in~. at iJ = and T = 0. At an arbitrary H 
the gap consists of two terms 



Vi = (l/4iV) ^(A/i + M2 - Mi_3 - M2-s)atatasai- 
V2 = (l/4iV)^(Mi - Ji)[2(a+a)i(a+a)_isin2a (10) 
— (a^a^)i(ai + a^Y) cos^ a], 

where 1 = qi etc. At small momenta we have in parentheses — 2A(3 • 4)/5' and Ak"^ / S for Vi and V2, respectively. 

Now one has to consider the Vi interaction and the second part of V2 interaction together. They give the principal 
contribution to the gap A^„j 

^lnt~ ^ 2.(A7q-Jq)- ^ 2.—, (11) 

where in r.h.s we take into account that ^ Jq = ^ Jci = ^ and according to Eq.(4) (k • c) = Ak"^ / [SDq). This 
contribution is T- independent. 

The contribution of the first term in V2 is more complicated one. Its T-independent part is proportional to (fca)^ 
and may be neglected [for MnSi {kaY ~ 0.03]. The T-dependent contribution consists of two parts. The excitations 
with q ^ k have a quadratic dispersion eq w Aq"^-- and at T ^ Ak"^ they are responsible for the first part, which 
has the form A| ^ = (l/2)(Afc2)2C(3/2)(fca)3[T/(27rAA:2)]3/2 sin^ acos^ a. In spite of small factor {kaY it may be 
important at sufficiently high T. 

The second part is not so trivial. According toi^ (cf. alsoii) at g < fc and H < He the spin-wave spectrum 
becomes strongly anisotropic due to umklapp processes connecting excitations with q and q ± k. In zero field taking 
into account the gap we have 

eq = [A\k\l + 3qi/8) + A'Y^'. (12) 

As the field increases this anisotropy becomes weaker since the term in the expression for eq appears to be proportional 
to q^. In the ferromagnetic state {H > He) this anisotropy vanishes. In a very weak field we can use Eq.(12) and 
then the second part is given by 

^2.2 = (l/7r)Afc2rv/3721n(AfcVA)sin2acos^a. (13) 

This expression diverges if A — > 0. Similar divergences were discussed in^^. However, we can neglect this contribution 
due to the small factor (kaY and decreasing of logarithm when H increases. 
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The three-point interaction is given by V3 = V_ cos a + V-^ sin a cos a where V± = {2S/Ny^'^J2Ci^\a+a)^^{a^^± 
a+) and 

= [Jt^ + ^q(q • c)] ^ (A/5)(q . k){qaf 
= [( - 4+))/2 - i?q(k . c)] ^ -AfcV(25), 

where the r.h.s. expressions are derived from Eq.(4) at g <C 1/a. The corresponding contributions to were 
evaluated as in2^. The T-independent term is proportional to (fca)^ and the second one is equal to — 2A2,2/3. Both 
of them are small and can be neglected. Finally we have 

A' = Aj„, + Al^,. (15) 

The field dependence of the ratio A^(iJ)/A^(0) for the three directions (111), (100) and (110) are shown in Fig. 3 
with He — 0.565r, A(0) = 0.1 T in agreement with the experiment. The cubic anisotropy A was chosen to be 
equal to —0.05 T. In case of H || [100] A^ remains positive and the spin-wave spectrum is stable at all H. In this 
case at H > Hq along with ferromagnetic spin configuration the spin-wave components of the lattice spins has to 
remain rotating [see Eq.(l)] as was discussed ini^. For (111) direction A^ is negative at H > Hi ~ 0.72Hc and the 
spin-wave spectrum becomes unstable. Hence the long-range helical order demolishes and the corresponding Bragg 
peaks disappear and the scattering has to be spread around them. Such decrease of the intensity along [111] is shown 
in Fig.l in qualitative agreement with the theory. However, more detailed measurements are needed. A similar 
instability has to be along [110] also contrary to the expectations for H || [100]. 

The width of this scattering may be estimated from the condition eq = 0, as for larger q the spin- wave spectrum 
can not feel the disorder. Near Hi we have ^ (Aqk)'^ + A^ and the inverse correlation length of the disorder 
K = k{\A\/Ak^) < k. Close to He - Aq^is and k = fc(|A|/Afc2)i/2. 

This disordered state have a strong chirality, which is demonstrated by a constant polarization of the scattered 
neutrons in the whole field range below He (see inset in Fig.l). This polarization is determined as P = (/_ — 
+ 1+) and according to the general theory the ratio P/Pq, where Pq is the initial neutron polarization, is 
the chirality at given q. A strong drop of the polarization at He is a signature of the first order transition to the 
uniform ferromagnetic state with weak chiral fluctuations. 

In conclusion, we analyzed thoroughly the field behavior of the spin wave gap in the spin-wave spectrum of cubic 
helimagnets. It is shown that if the field is applied parallel to both [111] and [110] directions a partially disordered 
state has to take place at Hi < He- We demonstrated that this state appears when the square of the spin- wave gap 
becomes negative and the spin-wave spectrum unstable. We presented the first observation in MnSi of this partially 
disordered chiral state in the magnetic field. 

The work is supported in part by the RFBR (projects No 05-02-19889, 06-02-16702 and 07-02-01318) and the 
Russian State Programs "Quantum Macrophysics" and "Strongly correlated electrons in Semiconductors, Metals, 
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FIG. 1: The intensity of the Bragg reflection in MnSi at T = 15 K as function of the field at H || [111]. The full line is the 
theoretical prediction (see the text). Inset: the spin chirality as a function of the field measured by polarized neutrons(see 
text). 




FIG. 2: Hartree-Fock (a) and three-point diagrams for the spin- wave gap (b). 
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FIG. 3: The magnetic field dependence of the ratio A^(_H')/A^(0). Parameters are given in the text. 



